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ABSTRACT 


Estimates of the average surface temperature on Mars are derived 
from radiative equilibrium considerations. A minimum possible surface 
temperature is estimated by computing the radiative equilibrium tempera- 
ture that the Martian surface would have if the planet had no atmosphere. 
An estimate of the maximum possible value of the average surface tempera- 
ture is obtained by computing the surface temperature that would result 
from a maximum greenhouse model. The computations indicate that the 
average surface temperature is in the range 219K to 233K. Comparisons of 

the theoretical computations with indications of surface temperature 

* 

obtained from thermal emission observations are found to be in reason- 
able agreement. 
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THEORETICAL ESTIMATES OF THE AVERAGE 
SURFACE TEMPERATURE ON MARS 


I. INTRODUCTION 

Of extreme interest to meteorologists and astrophysicists who are 
involved in planetary atmospheres research is the average surface tem- 
perature of the planets o Planetary surface temperature can be derived 
from both experimental observation and theoretical reasoning* The ob- 
servational determinations 3 which work well in the case of a rather 
tenuous atmosphere such as Mars** are based upon measurements of the 
planet's emission of infrared radiation; the theoretical determinations 
are based largely upon some sort of radiative equilibrium considerations* 
In the present study we derive estimates* based upon radiative equili- 
brium considerations* of the minimum and maximum possible value of the 
average temperature of the Martian surface* These estimates are then 
compared to some of the thermal emission observations* 
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II. THEORY 


The average surface temperature of a planet which has an atmosphere 
that is transparent to long wave radiation can be computed by equating 
the energy received from the sun to the energy emitted by the planets* 
surface. The following formula for the temperature of such a plant 
can be easily derived: 




If one now introduces into such a planetary atmosphere^ gases that 
absorb long wave radiation^ a "green-house 11 effect is created and the 
average surface temperature will increase above that given by Equation 
(l)o Thus,, for an atmosphere that is not completely transparent to 
long wave radiation Equation (1) can be used to derive an estimate of 
the lowest possible value of the average planetary surface temperature,. 
For the planet Mars the following numerical values can be substituted 
in Equation (1) c 
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A = 0.15 


M- = 0.44 

n -2 , "I 

s.Co = 2,0 cal crn mm 

-11 -2 ”4 -1 

cr = 8.13 x 10 cal cm deg min 

With these substitutions the lowest possible value for the average 

planetary surface temperature for Mars turns out to be 219K« 

The problem of estimating the maximum possible value of the average 
Martian surface temperature is more complex. We wish to maximize the 
green-house effect,, within reason, in order to obtain a realistic esti- 
mate of the maximum surface temperature. The green -house effect depends 
essentially upon the amount and vertical distribution of absorbing 
gases 3 and on the vertical distribution of temperature. The greater 
the amount of absorbing gases - the more opaque the atmosphere is - 
the greater will be the green-house effect,, and thus the higher will be 
the surface temperature. In our model (discussed later) we attempt to 
make realistic estimates of the maximum amounts of carbon dioxide., water 
vapor and ozone in the Martian atmosphere, and their distribution with 
height. Once the maximum amount of absorbing gas is fixed we have only 
to choose a realistic temperature distribution with height,, which would 
maximize the surface temperature. For simplicity,, and for lack of know- 
ledge,, we shall assume that the temperature variation with height is 
linear. This is probably a good approximation in the troposphere of 
Mars 3 as it is for the troposphere of the earth P and it is the tropo- 
sphere rather than the upper atmosphere that contributes most to the 
radiation fluxes. The question we must now answer is, "What choice of 
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constant lapse rate will produce a maximum green-house effect and thus 
a maximum surface temperature?" We can answer this question by going 
through the following qualitative reasoning* 

The upward flux of long wave radiation at the top of Martian 
atmosphere must balance the incoming solar radiation that is not re- 
flected back to space. This can be written as 

I 0 (l-A) = W (2) 

where is the average incoming flux of radiation at the top of 

the Martian atmosphere^ 

A is the planetary albedo^ and 

W is the upward flux of long wave radiation at the top of 
the atmosphere. 

The upward flux W is composed of radiation coming directly from the 
planet's surface and of radiation emanating from the planet* s atmosphere. 
If the atmosphere is isothermal $ the radiation originating in the atmos- 
phere will be at the same temperature as the radiation originating from 
the surface, Thus 3 the upward flux will remain constant with height 
and have a value equal to the black body flux emitted by the surface <, 

The surface temperature of an isothermal atmosphere would then be the 
same as that computed from Equation (1)„ i.e, 219K* If the temperature 
decreases with height, however^ the radiation, emitted by the atmosphere 
would be at a lower temperature than that emitted by the ground* In 
order that W now balance the incoming radiation,, the surface temperature 
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must be higher than 219K U The greater the lapse rate the greater is 
the surface temperature in order for balance to occur at the top of the 
atmosphere* In a planetary atmosphere the limiting lapse rate is the 
adiabatic lapse rate and we choose this as the lapse rate that will 
maximize the green-house effect and hence the surface temperature* 


With a given distribution of absorbers and temperature lapse rate, 
the outgoing radiation, W, is a function only of the surface temperature, 
and can be written as 


T 

p s 


T 

p s 


T 

p s 


W = cr T 


R(C0 2 ) dT 


R(H 2 0) dT 


R(0 3 ) dT (3) 


where R 


*2 

v 

d T 


(1 « t ) dv 


(4) 


and T^ is surface temperature, 

T^ is the temperature at the top of the atmosphere 
B is the black body flux, 
t ^ is the flux transmissivity 2 and 
v is frequency* 

R is a function of path lengthy u ? and temperature, T, and the particu- 
lar gas in question* Tables of R for carbon dioxide,, water vapor and 
ozone have been presented by Elsasser"^ and expanded by Hales,, et ah 
R for carbon dioxide covers the 15m band, R for ^0 covers the 6*3jJ. band, 
the window region,-, and the rotational band, and R for ozone covers the 
9o6p band. The integrals in Equation (3) can be evaluated numerically, 
given the temperature and absorber distribution in the atmosphere* It 
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should be noted that both water vapor and carbon dioxide absorb radia- 
tion in the 1 5\i band* However-, since we are attempting to maximize the 
green-house effect we shall treat these absorptions independently and 
make no correction for the overlap * 


Equation (2) can now be written as: 


T 

n s 


T 

p s 


T 

p s 


i„d - a> + 


R(CO ) dT + / R(H 0) dT + 


R(0 ) dT = cr T 
j s 


(5) 


To solve this equation for the surface temperature, T 3 one can 

s 

assume a value for and perform the integrations on the left side of 
the equation,. Upon adding the incoming radiation to the evaluated in- 
tegrals, we can compute a second approximation of from Equation (5)« 
This new value is then used to recalculate the integrals and obtain a 

third approximation, and so on, until a value of T is obtained that 

s 

balances Equation (5)* The surface temperature obtained from this model 
represents an estimate of the maximum possible value of the average 
surface temperature on Mars a 
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Ill, GREEN- HOUSE MODEL AND COMPUTATIONAL TECHNIQUES 


The green-house model that we shall adopt is based upon reason- 
able estimates of the maximum amounts of carbon dioxide , water vapor and 
ozone in the atmosphere of Mars* Of these three constituents, carbon 
dioxide is the only one that has definitely been detected in the Mirtian 
atmosphere, Grandjean and Goody in a re-analysis of Kuiper's^^ 
near infrared measurements 3 derive a carbon dioxide content of about 
2 % for the commonly accepted surface pressure of 85mb* We shall simply 
double this value for our estimate of maximum carbon dioxide content,. 

We shall further assume that all three gases are uniformly mixed with 
height so that the path lengths are a function only of pressure, and 
that a linear pressure correction can be applied directly to the path 
lengths as suggested by Elsasser* With these assumptions and a car- 

bon dioxide content of 4% the carbon dioxide path length for any pres- 
sure thickness is 


u co 2 - 86 • 


( — 

yooo/ 


( 6 ) 


where A U is the pressure corrected path length in a layer with mean 

_ 2 

pressure p mb and thickness A p t , 


Although water vapor has not been detected spectroscopically in 
the Martian atmosphere^ the aqueous nature of the polar caps indicates 
that there must be some water vapor present, Dunham^ ^ on the basis of 
spectroscopic observations^ derived an upper limit of 0o0015 times the 
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terrestrial amount over Mt. Wilson on an average clear night. This is 

of the order of 10" 3 cm of precipitable water with the terrestrial amount 

(6 ) 

of about Oo 7 cm prevailing during Dunham’s observations . DeVaucouleurs 
states that the amount of precipitable water is very likely much less 

o -2 

than 4 x 10 cm and probably less than 1 x 10 cm. We take as our esti- 

-2 

mate of the maximum water vapor content the value 10 cm of precipitable 
water,. Assuming a constant mixing ratio with height* the pressure cor- 
rected water vapor path length for any pressure thickness can be written 
as 


v 
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(7) 


Ozone also has not been detected in the Martian atmosphere . An 

(4) 

upper limit of 0.05 cm has been given by Kuiper*> who based his esti- 

mate on the failure of attempts to detect ozone by means of absorption 
spectra* Another estimate of the maximum amount of ozone has been de- 
rived by Mar mo 2 e t al , who computed the total amount of photochemi- 
cally produced ozone for an atmosphere that had a maximum possible 
amount of oxygen* The maximum amount of oxygen is about 2*5 m (Dunham 
and for this oxygen content Marmo ? et al computed a value of about 0*15 
cm for the total amount of ozone* Using this latter value as our esti- 
mate of the maximum amount of ozone on Mars we can express the pressure 
corrected ozone path lengths as 


) 


= 1.8 x 10 


(— E— ) 

Viooo / 


A p 


( 8 ) 
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As discussed previously the adiabatic lapse rate 5 approximately 
3 B 7 deg/km 5 is to be used in the radiation computations* The surface 
pressure is taken as 83 mb. In order to evaluate numerically the inte- 
grals appearing in Equation (5) P the atmosphere is divided into nine 
layers of thickness 10 mb^ except for the topmost layer , whose thickness 
is 5 mbo For example the integral for carbon dioxide is written as 

X \ Clog u i3 T £ ) a T. , 
i=l 

t h 

where is the R value of the i layer, which is a function of u^ 3 

the pressure corrected carbon dioxide path length from top of the atmos- 
phere to the middle of the layer and T^ 3 is the average temperature of 
th 

the i layer 5 and A Th is the temperature difference between the bottom 
t h 

and top of the i layer * 

Mean temperatures and path lengths for each layer are used to ob- 
tain values of the R function from Elsasser * s tables-, and the com- 
putation is similar in all respects to that suggested by Elsasser except 
that numerical rather than graphical integration is used* In our model 
the temperature at the top of the atmosphere is zero K and the tempera- 
tures in the upper part of the atmosphere are below 193K* Since Elsasser 7 s 
tables of R cover only the range 193 ^ T ^ 3 13 3 they were extrapolated 
from 193K to a temperature of zero K. for selected path lengths*, in our 
computations* The extrapolation to zero degrees K was performed with 

the knowledge that R at zero degrees is equal to zero (see Equation (4))* 

n 193 

As a check on the extrapolation., integrals of the form / R dT at 

4 > 
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constant path length were evaluated graphically and compared to similar 
integrals given by Elsasser. ^ ^ Although there may still be some error 
in the values of R used at temperatures less than 193K ? these will not 
significantly affect the computations since the contributions of the 
upper layers to the outgoing radiation are much less than the contribu- 
tions from the lower layers u Graphs showing the extrapolated R curves 
are shown in Figures 1, 2 ? and 3, 
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R(C0 2 ) - (cal -cm 2 day ') 



O too 200 


TEMPERATURE (°K) 


Figure 1. Extrapolated values of R(CO^). 
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R(H 2 0) - (ca/-cm~ 2 day~') 


0(6 A3 - 5 20 



TEMPERATURE (°K) 


Figure 2. Extrapolated values of RCH^O). 
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R( OjJ - (cal cm 2 day ') 


0.08 


Ql G A3- 530 



Figure 3. Extrapolated values of R(O^). 
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IV. RESULTS 


_2 -1 

For a solar constant of 2.0 cal cm min and a dilution factor of 

0.44; the average incoming radiation at the top of the Martian atmos- 

-2 -I 

phere^ I 3 is 317 cal cm day . With the commonly accepted value of 
0.15 for the Martian planetary albedo the average amount of solar radia^ 
tion remaining after reflection is 

I (1~A) = 269 cal cm ^day 1 .. 

Using this value and the maximum green-house model discussed above 9 we 
obtain from Equation (5) a value of equal to 233K„ This temperature 
is our estimate of the maximum possible value of the average surface 
temperature on Mars« Comparing 233K with our previous estimate of 219K 
for the minimum possible value 5 we see that even a maximum green-house 
effect would raise the surface temperature only about 14 deg K, This 
maximum green^house effect can be compared to the average green-house 
effect for the earth. The temperature of the earth c s surface^ if it 
had an atmosphere transparent to long wave radiation,, can be computed 
from Equation (l)* With an albedo of 0.35^ U - 1 and a solar constant 
of 2 o 0 cal cm ^min \ Equation (1) yields a temperature of 252K for 
the ear the The observed average temperature of the earth is about 288K P 
indicating that the average green-house effect for the earth is about 
36K 2 or about 2^ times as large as the maximum green-house effect for 
Mar s 


14 



Of the three gases contributing to the green-house effect, carbon 
dioxide is by far the most important as can be seen from the following 
values of the integrals appearing in Equation (5): 

T s 

I R(CO ) dT = 57 cal cm ^day ^ 

j T 2 

t 

T 



R(H 2 G) dT - 18 cal cm" 2 day 


RCO^) dT = 4 cal cm 2 day 


This is in contrast to the earth^s atmosphere where water vapor is the 
most important absorber of long-wave radiation 

Another measure of the effectiveness of the green~house is the 
ratio of outgoing long-wave radiation at the top of the atmosphere to 
the long-wave radiation emitted by the surface^ which may be termed the 
infrared transparency of the atmosphere,. The lower the infrared trans- 
parency,, the more effective is the green house*. For average conditions 

on earth the infrared transparency (based on radiation fluxes computed 

(9) 

by Landon * y ) is about 57%; for our maximum green house on Mars^ the 
infrared transparency is 77%* 

Based upon our theoretical estimates the true average surface 
temperature of Mars should lie between 21.9K and 233 K* In order to fix 
the average temperature more accurately^ it would be necessary to know 
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the actual amounts and distribution of absorbing gases 9 and the actual 


vertical temperature distribution. A sample calculation indicated that 
the use of the estimated average carbon dioxide concentration - 2% by 
volume - rather than a maximum concentration - 47o by volume - did not 
affect the outgoing radiation too much* Also the vertical temperature 
distribution we used cannot be too far from the actual temperature dis- 
tribution* We therefore suggest that the actual mean temperature is 
closer to 233K than to 219K e In the next section we compare our esti- 
mate with some of the observations of Martian temperatures. 
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Vc COMPARISON WITH OBSERVATIONS 


The observations of surface temperature on Mars are generally based 
upon measurements of the planet 5 s thermal emission; these measurements 
are usually made during oppositions and therefore refer largely to 

/Q \ 

the sunlit side of the planet* Gifford uss analyzed systematic** 

ally large numbers of such observations and has derived average latitudinal 
surface temperature distributions on the noon meridian for each of the 
four seasons for the latitude range 60°N to 80°S* On the basis of con- 
tinuity and solar insolation considerations we have extrapolated these 
curves to the poles; they are shown in Figure 4$ with extrapolations 
indicated by dashed lines* In order to compute the average planetary 
surface temperature from these observations we must perform the follow- 
ing operations: 

(1) Average the seasonal curves to obtain an average 

annual curves 

(2) Subtract from these noon^time temperatures the am- 
plitude of the average diurnal variation of temperature as a function 
of latitude in order to obtain a mean daily temperature; 

(3) Compute the average planetary surface temperature 
by neigMii'g the mean temperature at: each latitude belt according to 
the area of the latitude belt* 

When the first of these operations is performed the average annual 
noon curve shown in Figure 5 results* The correction of this curve for 
diurnal variation is important because such variations are apparently 
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Figure 5. Average annual surface temperature as a function 
of latitude ; (a) noon temperatures (b) mean daily 



appreciable on Mars u There have been a number of measurements;, mostly 
in tropical regions, of the variations of Martian surface temperature 
during the day w Gif ford 5 s V ' analysis of 214 daytime surface tempera- 
tures near the Martian equator indicates a diurnal range of about 70C, 
which is equivalent to an amplitude of 35C r DeVaucouleur s in his book 
on Mar s , suggests an average amplitude of 25C in the equatorial regions « 
In the recent observations of Sinton and Strong however, the indi- 

cated amplitude is about 45C or higher „ Since the analyses of Gifford 
and DeVaucouleur s are based upon more data we shall rely on their esti- 
mates for the amplitude of the diurnal variation at the equator; we 
therefore adopt an amplitude of 30C in the equatorial regions of Mars, 
Assuming that the diurnal variation depends only on insolation, we can 
approximate the latitudinal variation of the amplitude by a simple 
cosine function 

a cp = ^ C0S ^ (9) 

where a^ is the amplitude at any latitude cp* When these amplitudes are 
subtracted from the noon curve, the average annual curve shown in 
Figure 5 is obtained., The curve indicates an average equatorial tempera- 
ture close to 2.40K and average polar temperature of between 200K and 
210K„ The average planetary surface temperature can be obtained from 
a numerical integration of this curve from pole to pole, weighting the 
temperatures according to surface area as follows? 

_ [ r n/2 - 

T = j / T (cp) cos cp d cp (10) 

- rt /2 
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where T is the average planetary surface temperature 0 and T (cp) is the 
average surface temperature at latitude cp„ Upon performing this inte- 
gration we obtain a value of 233K for the average Martian surface tempera- 
tures, Within the scope of the uncertainties involved in obtaining this 
value, it can be considered to be in reasonable agreement with the com- 
puted theoretical range 219K to 233K* 
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VI. SUMMARY 


From theoretical computations based upon radiative equilibrium 
considerations, the average surface temperature of the planet Mars is 
found to be in the range of 219K to 233K, 

This is in reasonable agreement with the thermal emission observa~ 
tions which suggest a mean temperature of about 233K U The computations 
also indicate that the maximum green-house effect on Mars is about 40% 
the .average green-house effect in the earth's atmosphere, and that 
carbon dioxide is the most important contributor to the Martian green 
house* 


22 



REFERENCES 


Elsasser. W„ M„, Atmospheric Radiation Tables, Amer. Meteor. Soc „ 
Meteorological Monographs .. 4, 23., 43 (I960) „ 

Hales, J, V., Williams, T. L. , Henderson, D. , "Calculation of 
Infrared Radiative flux Emission of the Earth Plus Atmosphere at 
Various Levels High /' t'-e the Earth". Final Report Contract AF 
19 (604) -2418. (1960). 

Gran.djean, J . , and Goody, R. M. , "The Concentration of Carbon 
Dioxide in the Atmosphere of Mars", Ap. J., 121, 548 (1955). 

Kuiper, G. P.., T he Atmospheres of the Earth and Planets . Univ. of 
Chicago Press, 434 (1952). 

Dunham, T„ ."Spec troscopic Observations of the Planets at Mount 
Wilson", The Atm ospheres of the Earth and Planets, Ed. bv G„ P 
Kuiper, 288=305 (1952). 

DeVaucouleurs. G„, Physics of the Planet Mars. Faber and Faber 
Limited, London, 365 (1954) „ 

Marino, F. F, , Watanabe, K. , Dalgarno, A., Warneck, P., "Laboratory 
and Theoretical Studies in the Vacuum Ultraviolet for the Inves- 
tigation of Chemical Physics of Planetary Atmospheres", Third 
Quarterly Progress Report, Contract No. NASw-124, GCA, Bedford, 
Mass., 70 (1960). 

Gifford, F„, "The Surface Temperature Climate of Mars", An. J 
123 s 154-161 (1956). ^ 5 

Landon, J., , "A Study of the Atmospheric Heat Balance", Final 
Report Contract No. AF 19(122)-165, AD 117227, Dept, of Meteorology 
and Ocenography, N.Y.U., 99 (1957). 

Sinton, W. M„, and Strong, J., "Radiometric Observations of Mars" 
AP* 131, 459-469 (1960). 



